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I.  SUMMARY 


This  final  technical  report  describes  research  work  ‘'Theoretical  and 
Numerical  Study  of  Anomalous  Turbulent  Transport  in  Plasmas"  performed 
under  AFOSR  contract  86-0156  during  the  period  of  June  15,  1986  to  December 
14,  1990. 

Theoretical  and  computational  analysis  of  particle  transport  in  a  plasma 
turbulence  was  carried  out.  Particle  behavior  was  investigated  by  test  particle 
simulations  in  plasma  turbulences,  including  drift  wave  turbulence,  ion  acoustic 
turbulence  and  Langmuir  turbulence.  A  semiclassical  approach  to  analyze  plasma 
turbulence  was  taken  to  understand  particle  transport. 

The  report  includes  a  brief  description  of  research  accomplishments  under 
the  support  of  this  grant  in  Sec.  II.  Section  m  describes  the  undergraduate 
fellowship  program  associated  with  this  grant.  Section  IV  lists  conference 
presentations,  while  Sec.  V  lists  papers  published  in  professional  journals.  Section 
VI  lists  the  titles  of  master  theses  and  PhD  dissertations  supported  by  this  grant. 
The  report  concludes  in  Sec.  VII  by  listing  staffs  involved  in  this  contract. 
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II.  RESEARCH  ACCOMPLISHMENTS 


Drift  wave  turbulence  — —  Anomalous  transport  of  electrons  in  the  low 
frequency  electrostatic  turbulence  (drift  wave  turbulence)  has  been  studied 
extensively.  Electrons  were  known  to  diffuse  across  a  magnetic  field  by  E  x  B 
drift,  where  E  is  a  turbulent  electric  field  fluctuation  and  B  is  a  background 
magnetic  field.  We  consider  the  transport  of  electrons  in  the  presence  of  shear 
structure  in  the  magnetic  field.  We  have  formulated  the  transport  of  electrons  in 
a  moderately  strong  turbulence  based  on  the  drift  kinetic  equation.  Electrons  are 
in  resonance  with  waves  at  a  mode  rational  surface,  which  is  determined  by  the 
rational  numbers  associated  with  Fourier  modes  of  drift  waves.  We  have  shown 
that  the  resonance  region  broadens  as  the  turbulence  becomes  stronger  and  that 
the  broadening  develops  with  time.  Our  new  formulation  of  resonance 
broadening  explains  our  observation  of  particle  transport  in  the  test  particle 
numerical  experiment.  The  resulting  transport  coefficient  can  describe  nonlinear 
aspect  of  wave-particle  interaction  in  the  drift  wave  turbulence.  This  principal 
investigator  was  invited  to  give  a  talk  on  this  subject  at  the  1989  IEEE 
International  Conference  on  Plasma  Science  at  Buffalo,  New  York  in  May  1989. 
The  full  paper  is  published  in  Physics  of  Fluids  (1990). 

Langmuir  Turbulence  —  Time  dependence  of  velocity-diffusion 
coefficient  of  plasma  particles  in  Langmuir  turbulence  has  been  studied  based  on 
the  generalized  Langevin  equation  and  the  derived  Fokker-Planck-type  equation. 
The  time  dependent  nature  of  the  diffusion  coefficient  is  caused  by  the  retarded 
effect  of  the  turbulent  collisions,  which  is  responsible  for  the  non-Markovian 
stochastic  process  in  a  Gaussian  (or  non-Gaussian)  system.  The  projection 
operator  method  has  been  used  to  derive  the  evolution  equation  known  as  the 
generalized  Langevin  equation  for  the  particle  velocity  in  an  electrostatic  plasma 
turbulence.  Upon  completion  of  hundreds  of  realizations  in  a  test  particle 
simulation,  we  obtained  numerical  results  in  agreement  with  analytical  results. 
Our  results  were  reported  by  a  series  of  papers  at  IEEE  International  Conference 
on  Plasma  Science  and  at  the  Annual  Meeting  of  Division  of  Plasma  Physics  of 
the  American  Physical  Society. 

Ion  Acoustic  Turbulence  —  Nonresonant  interaction  of  plasma  particles 
and  plasma  waves  has  been  studied  in  the  presence  of  ion  acoustic  turbulence.  We 
have  confirmed  by  test  particle  computer  experiments  that  particles  away  from 
the  resonance  layer  could  indeed  diffuse  as  was  predicted  by  our  analytical  work. 
The  results  were  reported  in  the  Annual  Meeting  of  Division  of  Plasma  Physics 
of  the  APS  (1987). 
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Relativistic  particle  acceleration  —  Relativistic  particle  acceleration  by 
plasma  turbulence  has  been  studied  in  a  quasilinear  regime.  Relativistic  charged 
particles  can  be  accelerated  by  taking  energy  from  plasma  turbulence.  We  have 
formulated  the  problem  classically  and  quantum  mechanically.  The  theory 
describes  the  interaction  between  charged  particles  and  plasmons.  The  self-field 
created  by  the  particle  itself  is  also  considered.  Our  equation  reduces  to  the 
Balescu-Lenard  equation  in  the  classical  limit  This  result  was  presented  at  the 
APS  Topical  Conference  on  Plasma  Astrophysics  at  Santa  Fe  in  September,  1988. 

Turbulence  in  toroidal  geometry  —  Analytical  study  of  turbulent 
fluctuations  in  toroidal  geometry  has  been  carried  out  in  collaboration  with 
Professor  A.  Hirose  of  University  of  Saskatchewan,  Canada.  In  a  collisionless 
plasma  we  have  generalized  the  ideal  MHD  mode  equation  to  describe  drift- 
Alfven  mode  and  examined  numerically  the  stability  problem  (published  in 
Nuclear  Fusion  (1987)).  In  a  collisional  plasma,  the  resistivity  could  be 
responsible  for  the  instability.  We  have  examined  the  resistive  MHD  balloning 
mode  by  solving  the  second  order  mode  equation  (published  in  the  Canadian 
Journal  of  Physics  (1988)  and  Nuclear  Fusion  (1989)).  The  mode  equation  is 
extended  to  incorporate  the  plasma  inhomogeneity  like  ion-temperature-gradient. 
The  numerical  study  of  the  ion-temperature-gradient  instability  is  published  in 
Nuclear  Fusion  (1987). 

Semiclassical  approach  to  plasma  turbulence  —  We  have  applied  the 
method  of  quantum  electrodynamics  to  the  analysis  of  plasma  turbulence.  The 
wave-particle  interaction  is  viewed  as  a  scattering  process  between  plasma 
particles  and  quasiparticles(plasmons,  phonons).  The  particle  transport  in  plasma 
turbulence  is  then  closely  related  to  the  transition  probability  of  absorption  and 
emission  processes.  Successful  application  of  the  method  resulted  in  clarification 
of  the  nonresonant  wave-particle  interaction  in  plasma  turbulence  (published  in 
Physical  Review  A  (1987)).  On  the  other  hand,  such  a  viewpoint  is  applied  to  the 
interaction  between  photons  and  plasmons,  and  a  novel  scheme  of  photon 
acceleration  (or  frequency  up-shifting  of  electromagnetic  wave)  is  proposed 
(accepted  for  publication  in  Physical  Review  A). 
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III.  UNDERGRADUATE  RESEARCH  FELLOWSHIP  PROGRAM 


The  undergraduate  research  fellowship  program  was  initiated  by  Dr.  Robert 
Barker  of  AFOSR  and  our  plasma  laboratory  was  involved  in  the  program 
during  the  period  of  the  contract.  The  selected  undergraduate  students  were 
involved  in  the  related  research.  They  participated  in  a  weekly  plasma  seminar 
with  graduate  students  and  presented  their  assigned  topics  at  the  seminar.  The 
names  of  undergraduate  students  and  the  titles  of  their  final  reports  are  listed 
here. 

From  Mav  16  to  August  31.  1987  : 

1.  Theodore  Christopher  Grabowski 

"Anomalous  Particle  Transport  in  Drift  Wave  Turbulence  in  a  Plasma" 

2.  Christopher  James  Koath 

"Study  of  Diffusion  and  Friction  in  an  Electrostatic  Turbulent  Plasma" 

3.  David  Edward  Stahlke 

"Ion  Acoustic  Plasma  Turbulence  and  Data  Processing" 


From  June  15.  1988  to  Mav  31.  1989  : 

1.  Hok  Tung  Wong 

"Application  of  MAGIC  Code  to  die  Study  of  Plasma  Turbulence" 

2.  David  Stahlke 

"A  Computational  Study  of  Various  Two-Stream  Instabilities" 

3.  Mark  Thompson 

"A  Numerical  Experiment  of  Relativistic  Particle  Diffusion  in  a  Plasma 
Turbulence" 
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rv.  LISTING  OF  CONFERENCE  PRESENTATIONS 

Conference  presentations  of  the  research  results  obtained  under  the 
contract  AFOSR  86-0156  are  listed  here. 

1.  Nonresonant  Wave-Particle  Interaction  from  the  Quantum  Mechanical 

Viewpoint,  O.  Ishihara,  Bull.  Am.  Phys.  Soc.  31, 1560  (1986). 

2.  Numerical  Study  of  Tuibulence,  E.  Robinson  and  O.  Ishihara,  Bull.  Am. 

Phys.  Soc.  31, 1497  (1986). 

3.  Real  Space  Diffusion  in  Drift  Wave  Turbulence,  E.  Robinson  and  O. 

Ishihara,  IEEE  Int.  Conf.  Plasma  Sci.,  paper  2X5  (Arlington,  VA  1987). 

4.  Analytical  and  Numerical  Study  of  Anomalous  Friction  in  a  Plasma 

Turbulence,  W.  Ho  and  O.  Ishihara,  IEEE  Int.  Conf.  Plasma  Sci.,  paper 
5T3  (Arlington,  VA  1987). 

5.  Resistive  Modes  in  Tokamaks,  A.  Hirose,  T.L.  Kroeker,  and  O.  Ishihara, 

Bull.  Am.  Phys.  Soc.  32, 1773  (1987). 

6.  Study  of  Velocity  Diffusion  of  Test  Particles  in  Ion  Acoustic  Turbulence, 

K.M.  Yuen,  O.  Ishihara,  and  A.  Hirose,  Bull.  Am.  Phys.  Soc.  32,  1859 
(1987). 

7.  Quasilinear  Effects  on  Turbulent  Bremsstrahlung,  O.  Ishihara  and  A. 

Hirose,  bull.  Am.  Phys.  Soc.  32,  1859  (i987). 

8.  Stochastic  Heating  of  Plasma  Particles  in  Ion  Acoustic  Turbulence,  K.M. 

Yuen  and  O.  Ishihara,  Spring  Meeting  of  the  Texas  Section  of  the 
American  Physical  Society  (Austin,  TX  March,  1988). 

9.  Plasma  Particle  Diffusion  in  Drift  Wave  Turbulence,  C.  Grabowski  and  O. 

Ishihara,  Spring  Meeting  of  the  Texas  Section  of  the  Am.  Phys.  Soc. 
(Austin,  TX  March,  1988). 

10.  Stochastic  Diffusion  in  the  Absence  of  Resonances,  O.  Ishihara,  W.  Ho, 

K.M.  Yuen,  and  A.  Hirose,  IEEE  Int.  Conf.  Plasma  Sci.  (Seattle,  WA 
1988). 

11.  Resistive  MHD  Ballooning  Mode  in  Tokamaks,  A.  Hirose  and  O.  Ishihara, 

IEEE  Int.  Conf.  Plasma  Sci.  (Seattle,  WA  1988). 

12.  Relativistic  Quasilinear  Particle  Acceleration  in  Plasma  Turbulence,  O. 

Ishihara,  APS  Topical  Conf.  on  Plasma  Astrophysics  (Santa  Fe,  NM 
1988).  Bull.  Am.  Phys.  Soc.  34,  1287  (1989). 

13.  Saturation  of  Nonlinear  Drift  Mode  Induced  by  Toroidicity,  O.  Ishihara 

and  A.  Hirose,  Bull.  Am.  Phys.  Soc.  33,  1941  (1988). 

14.  Time-Dependent  Diffusion  Coefficient  for  Non-Markovian  Process  in 

Plasma  Turbulence,  H.  Xia  and  O.  Ishihara,  Bull.  Am.  Thys.  Soc.  33, 
2020  (1988). 
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15.  Study  of  Stochastic  Heating  Mechanism  of  Plasma  Particles  in  the  Presence 

of  Ion  Acoustic  Turbulence,  K.M.  Yuen  and  O.  Ishihara,  Fall  Meeting  of 
the  Texas  Section  of  the  American  Physical  Society  (Lubbock,  TX  1988). 
Bull.  Am.  Phys.  Soc.  34,  1521(1989). 

16.  Velocity  Diffusion  Coefficient  as  a  Function  of  Time,  H.  Xia  and  O. 

Ishihara,  Fall  Meeting  of  the  Texas  Section  of  the  Am.  Phys.  Soc. 
(Lubbock,TX  1988).  Bull.  Am.  Phys.  Soc.  34,  1521  (1989). 

17.  Particle  Diffusion  in  Turbulent  Fields:  Transition  from  Quasilinear  to 

Nonlinear  Stage,  O.  Ishihara,  IEEE  Int.  Conf.  Plasma  Sci.  (Buffalo, 
NY  1989).  (Invited  Talk). 

18.  Generalized  Langevin  Equation  for  Plasma  Turbulence,  H.  Xia  and  O. 

Ishihara,  IEEE  Int.  Conf.  Plasma  Sci.  (Buffalo,  NY  1989). 

19.  Study  of  Properties  of  the  Stochastic  Processes  in  Plasma  Turbulence,  H. 

Xia  and  O.  Ishihara,  Bull.  Am.  Phys.  Soc.  34,  1926  (1989). 

20.  Photon  Acceleration:  Three  Wave  Interaction,  O.  Ishihara,  IEEE  Int.  Conf. 

Plasma  Sci.  (Oakland,  CA  1990). 

21.  Effect  of  Turbulent  Collisions  on  Diffusion  in  Stationary  Plasma 

Turbulence,  H.  Xia  and  O.  Ishihara,  IEEE  Int  Conf.Pasma  Sci.  (Oakland, 
CA  1990). 

22.  Photon  Acceleration  by  the  Plasmon-Photon-Photon  Interaction,  O. 

Ishihara,  Bull.  Am.  Phys.  Soc.  35,  2125  (1990). 

23.  Analysis  of  the  Non-Markovian  Behavior  in  Velocity  Diffusion  Process, 

H.  Xia  and  O.  Ishihara,  Bull.  Am.  Phys.  Soc.  35,  1956  (1990). 

24.  Photon  Acceleration  by  Plasma  Turbulence,  O.  Ishihara,  IEEE  Int.  Conf. 

Plasma  Sci.  (Williamsburg,  VA  1991). 

25.  Characteristics  of  Distribution  Function  in  a  Time-Dependent  Velocity 

Diffusion  Process,  H.  Xia  and  O.  Ishihara,  IEEE  Int.  Conf.  Plasma  Sci. 
(Williamsburg,  VA  1991). 
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V.  LISTING  OF  PUBLICATIONS 


Publications  of  the  research  results  obtained  under  the  contract  AFOSR  86- 
0156  are  listed  here. 

1.  Nonresonant  Wave-Particle  Interaction  in  Semiclassical  Quasilinear 

Theory,  O.  Ishihara,  Phys.  Rev.  A35,  1219-1225  (1987). 

2.  Drift  Alfven  Eigenmode  in  Tokamaks,  A.  Hirose  and  O.  Ishihara,  Nuclear 

Fusion  27,  1231-1238  (1987). 

3.  Search  for  Ion  Temperature  Gradient  Driven,  Electrostatic  Hydrodynamic 

Instability  in  Tokamaks,  A.  Hirose  and  O.  Ishihara,  Nuclear  Fusion  27, 
1439-1451  (1987). 

4.  Resistive  Ballooning  Mode  in  Tokamaks,  A.  Hirose,  T.L.  Kroeker,  and  O. 

Ishihara,  Can.  J.  Phys.  66,  1069-1075  (1988). 

5.  Resistive  MHD  Ballooning  Mode  in  Tokamaks,  A.  Hirose  and  O.  Ishihara, 

Nuclear  Fusion  29, 795-803  (1989) 

6.  Resonance  Broadening  in  Drift  Wave  Turbulence,  O.  Ishihara,  C. 

Grabowski,  and  A.  Hirose,  Phys.  Fluids  B2,  270-279  (1990). 

7.  Photon-Plasmon-Photon  Interaction,  O.  Ishihara,  Phys.  Rev.  A  (in 

press). 
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VI.  LISTING  OF  THESES  AND  DISSERTATIONS 


Master  theses  and  doctoral  dissertations  supported  by  the  contract  AFOSR 
86-0156  are  listed  here. 


Maser  Theses 

1.  Edward  Kyle  Robinson,  "Real  Space  Diffusion  in  Drift  Wave 
Turbulence"  (1987). 

2.  Kin-Ming  Yuen,  "Velocity  Space  Diffusion  in  Ion  Acoustic  Turbulence" 

(1989). 

PhD  Dissertations 

1.  Wai  H.  Ho,  "Wave-Particle  Interaction  in  Plasma  Turbulence"  (final 

version  is  in  preparation). 

2.  Huajuan  Xia,  "The  Non-Markovian  Process  in  Velocity  Diffusion  in  the 

Stationary  Plasma  Turbulence"  (in  preparation). 
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VI.  STAFFS 


Faculty  Investigator: 

Professor  Osamu  Ishihara  of  Department  of  Electrical  Engineering  and 
Department  of  Physics  at  Texas  Tech  University  is  a  principal  investigator. 

Graduate  Research  Assistants: 

During  the  period  of  the  contract,  two  doctoral  students  and  two  master 
students  were  supported. 

Undergraduate  Research  Assistants: 

Five  undergraduate  students  were  supported  under  the  program  of 
undergraduate  research  fellowship. 


Accountants: 

Financial  record  was  kept  by  Ms.  Sandra  Branch  of  the  Department  of 
Electrical  Engineering.  She  was  partially  supported  by  this  contract. 
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Program  of  the  Twenty-Eighth  Annual  Meeting 
of  the  Division  of  Plasma  Physics 
Baltimore,  Maryland;  3—7  November  1986 
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partlcla-quaslpartlcla  Intar  action  contributes  to 
transition  probability.  Tba  roosnsnt  transition 
probability  (1).  know  as  Fsral's  golden  rula.  Is  now 
supplemented  by  a  noma  sonant  pact,  which  guarantees 
tba  conssrrstlon  of  energy  and  momentum  In  tha  systea 
aa  Is  predicted  by  tha  classical  quaslilnear  theory. 

*  (fork  supported  by  AfQSE  under  grant  AfOSK-M-OlSd. 
(1)  C.  C.  Harris,  in  Advances  in  -Plasms  physics,  ed.  _ 
by  A.  Sleoe  and  V.  1,  thonpoon.  Vol.  3.  p.lS7 
(Ullsy.  Haw  York,  19*9) . 
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Munarlcal  Study  of  Toot  Part  tela  Tranonort  In  Prlft  ! 
Wave  Turbulence. «  E.  loblnaon  and  0.  Xahlhara,  Texas  1 
Tech  University  —  In  drift  vave  turbulence  where  the 
uagnet lc  field  is  expressed  as  I  •  to(0,x/L^.l),  L,- 
ahear  length ,  a  real  apace  diffusion  coefficient  ban 
been  formulated*  and  Is  found  to  hare  slactllar  char¬ 
acteristics  as  the  tine  dependent  diffusion  coefficient 
derived  for  velocity  space  in  Lancralr  turbulence.2 
1  numerical  experiment  reveals  the  particle  diffusion 
In  real  space  In  the  x  direction.  Test  particles  are 
fallowed  In  their  trajectories  in  the  prescribed  ran¬ 
dom  electrostatic  fields,  4(2.0  ,^_exp[l(kjr 

♦  k^r  -w^t)]  ,  where  f  Is  tha  electric  potential 
of  Che  (n,n)  node.  A  preliminary  result  of  the  ob¬ 
served  diffusion  coefficient  demonstrates  the  tint 
dependent  structure  of  the  diffusion  coefficient  so 
predicted  by  the  analytical  study. 

*  Wotfc  support ad  by  AFOSt  under  grant  AFOSt-86-01$6 

1.  A.  Salat.  Z.  Haturforach.  Tell  A  38,  1189  (1983). 

2.  0.  lohlhara  and  A.  Rlrose,  Phya.  Fluids  28.  21S9 

(1583),  — 
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Real  Spaoe  Diffusion  In  Drift  Wave  Turbulence* 

E.  Robinson  and  O.  tstrihara 

Department  of  Bectrtcal  Engineering 
Texas  Tech  University.  Lubbock,  Texas  79409 

In  drift  wav*  turbulence  characterized  by  a  sheered  mag¬ 
netic  held  6  -  B0(ez  +  x/Lsey),  where  L,  is  the  shear 

length,  and  random  electrostatic  fields.  $(x,t)  -  X  $mn 
X  +  Ni*  where  mn  Is  the  electriepotentiat 

of  the  (nun)  mode,  a  numerical  experiment  has  been  carried 
out  to  examine  the  particle  dfflusion  process  in  teal  space.  Two 
cases  are  examined,  one  with  an  Infinite  number  of  n-modes 
and  the  other  with  a  single  n-mode.  100  m-modes  are  used  lor 
both  cases.  For  infinite  n-modes,  diffusion  coefficients  found 

are  in  agreement  with  the  value  0  -  *]£  (|m$>m|2)  as  was 

predicted1  for  a  low  turbulence  level,  while  lor  a  single 
n-mode  w e  observe  that  diffusion  greatly  depends  on  the  choice 
of  mandn. 

An  analytical  expression  for  the  space  diffusion  in  the 
x -dir action  is  given  by 

-  ^  m2Dt2(3s  -  2t) 

X  e  3  (1) 

where  X  is  the  initial  position  of  a  test  particle,  extribits 
analogous  characteristics  to  the  qussffinear  velocity  space 
diffusion  coefficient.2  For  low  turbulence  levels,  Eq.  (1) 
shows  quasilinear-ftke  behavior. 


I 

J  « 


dt  e,(mx  +  ")i 
S-T 


D(X)-nl(|m$J2)8(rnX  +  n).  (2) 

at  the  resonant  locations.  X  -  -fnfm).  From  further 
examination  of  the  single  n-mode  at  a  rather  high  turbulence 
level,  we  have  observed  deviation  of  the  diffusion  coefficient 
from  that  of  Eq.  (2),  le.  high  turbulence  levels  cause  diffusion 
in  the  space  surroundng  the  resonant  locations  as  Is  predteted 
by  Eq.  (1).  A  comparative  study  between  the  real  space 
diffusion  in  drift  wave  turbulence  and  the  velocity  space 
rtiffusJon  In  Langmuir  turbulsnoe,  c j  suip,,2  at  a  rather  strong 
turbulence  level  demonstrates  fits  fime  dependent  nature  of  the 
diffusion  coefficient  as  was  described  by  Eq.  (1). 


*  Work  supported  by  AFOSR  under  grant  AFOSR-06-O156. 
1.  A.  Salat  Z.  Naturforsch.  Teti  a  3L  1189  (1983). 

Z.  O.  IsNhara  and  A.  Hboee.  Phya.  Fluids  2L  2159  (1985). 
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rm  i :  f  i 1 , ,  v ,  f-’i  s.'Tttt.  tt?  r.rr 


Wai  Ho  and  Osamu  Iduban.  Texas  Tech  University,  Lubbock. 
TX  79409-4439.-  »  -The  naociotouiihmit^  effect  oo  pltsnui 


^'zrrrZ'VTi  rr 


oi BUgpcoc ocm if  gmoioq  iiunirr iciiiy  m nmjraciBy.  antsoeen 
known  that  the  frictional  effect  together  with  die  effect 

electric  Gelds,  Sbeditaonion  of  panicle  trajectories  can  be  given  by 

dS<0  mjstiO  df  .where  £v(t)  «1g^6(S(0/l  dtf. 

We  thus  define  the  friction  coefficient  as*F^»4)*  d<A^0>Ak. 

It  ii  straightforward  to  show  itm  ifae  ttcdoBCOcffidmli^wi  by 

Ffo)  «^-)y A^/ds  Jx  .^(Vw)  <exp(£>  4x(t-i)> 

,y0  /  .  dv  (I) 

»  “t* 

where  *<*>-{  /|f-?  (teflexpUfo*^*! 
and 

Alternatively.  we  can  write  *P(v,t)  as 

f(v.t)  * 4-fi|r * /dr/'ds f* ,JX  <^.0^(0. 0>POf.vAs) 

^  37  ,'r.  P) 

1 

2  I  r  /JAi  *| 

where  Pp^xs)  -  (2k<a3<-s)]  >P  “PI'  2<SV>J  * 

A  numerical  experiment  has  been  carried  out,  in  which  one 

djmMidflnij  ytvlnm  >WfA»  fields  Bt  modeled  fof  LJOpPUlf 

turbulence  and  are  described  by  ^ 

<E(x,t)E(0,0>  «  -  [***&  ljxl)  cosCkex  -  «pt)] 

where  «pb  the  plasma  frequency.  kj  is  the  center  wiwiwnhrr  and 
q  la  the  nonnaHxadnna  value  of  the  turbulent  potential  lathe 
limiting  case  of  weak  nnbulenoe  (a  <<0.01),  we  have  observed  a 
well-defined  friction  coefficient  as  given  by  the  Both  SS—fO 
of  Eq.(l); 

^  « ■*$>  C3) 

la  the  modetaae  turbulence  level  (p  ■>  OjOI),  the  friedan  coefficient  is 


i  W  '  i  ’•  I  f  . 


to  dhow  some  broadening  effect  as  lep^^ed  by  <Ua  term  fat 
Eq.  O). 

The  formulation  of  the  friction  coefficient  b  also  studkd  by  a 
renonnatod  method  with  Green's  function  lysed  oo  theVlasov 
equation/  With  the  help  of  diagtaninuuic  technique,  we  formulate 
die  mbulent  friction  by  taking  secular  terms  fat  die  Green’S  function 
which  agrees  with  Eq.  (1)  in  some  limiting  situation. 

*  Work  supported  by  AFOSR  under  grant  AFOSR464136. 


1. 0.  Udhara  and  A.  Hirose.  Phys.  Ruids^,  2159(19*5). 

2.  M.  Kono  and  Y.  H.  Ichikawa,  Prog.  Theor.  Phys.JJ.754<1973). 


Program  of  the  Twenty-Ninth  Annual  Meeting 
of  the  Division  of  Plasma  Physics 
San  Diego,  California;  2-6  November  1987  . 


JR II  lulttlvc  HodM  la  »-■ - ■  a,  mmom, 

T-l~  noro-pnlT.’  ofta^:  Q.  UCTOiS  Tach 
B«l».  —Ik  iwImIn  |  «1  rlffUM  -trlii  la  tnlsrrt- 
— — *nr  bava  Wm  WMMlMl  la  tana  of  the  fantralUad 
■oaaata*  fcee  froa  aay  order  lag  aneog  frafaaadae,  w, 

V  -  JV.  <«>*»  ♦  an  alMO  (the  wagaatle 

fclft  fraqmacy  la  the  «t  apace).  Cadi  gaaacal 

iaatloa  la  anew ary  for  twalsctvw  wodas  trrriTi  of 
atzaac  tallal  local lxatfoa  ahoat  a  rarlnaal  aarfaca  and 
cocraapoadlat  largo  Made  width  la  n-apaca.  LkOOl  » 

H  »  “*■  *>*  typical  tokaaak  M(a  paraaatm.  a  aaocr- 
lcal  oearch  for  anatahla,  alpaafaactloaa  of 

raaladoe  nodea  has  fallal.  Only  the  atahla  (scahlllxad 
by  shear  and  toroldldty)  drift  aoda  (predicted  analyti¬ 
cally)  has  been  found. 

•Sponsored  by  KSERC  (Canada)  and  USAF-OS* 

1A.  Hlroee.  Hucl.  Fusion  (la  press).  Also.  Com. 
riasaa  Rqra.  Controlled  Pus  loo  10,  229  (1987) 
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Study  of  Velocity  Diffusion  <of  Test  Particles  in  Ion  Acoustic 
Turbulence.*  kIm.  Yuen  and  O.  Ishihara,  texas  Tech  Hnhrw 
and  A.  Hirose,  Unix.  of  Saskatchewan, — tn  order  to  study 
velocity  diffusion  3  ions  tn  ion  acoustic  turbulence,  we  set 
the  test  particles  in  prescribed  turbulent  electric  fields.  The 
numerical  experiments  are  carried  out  with  particles  placed 
randomly  in  the  system,  or  'alternatively  we  followed  the 
test  particles  with  the  same  initial  position  lor  different 
realization  of  electric  fields.  The  ion  acoustic'  turbulence  is 
modeled  by  the  dispersion  relation  tMch  is  characterized  by 
either  nondispersive  or  dispersive  nature.  The' results  of  the 
numerical  experiments  indicate  that  in  both  nondispersive 
and  dispersive  cases,  ions '  which  do  not  satisfy  the  resonant 
condition  do  diffuse  In  velocity  space.  The  analytical  study- 
also  shows  a  possibility  of  diffusion  of  nonresonant  particles 
if  the  turbulent  spectrum  has  a  finite  correlation  time.' 
The  numerical  results  show  fair  agreement  with  theoretical 
predictions  provided  that  the  initial  velocities  are  set  within 
critical  values,  which  are  located  well  above  and  well  below 
the  phase  velocity  of  the  waves.  It  has  also  been  observed 
that  there  is  no  stochastic  heating  of  particles  if  their  ini- 
tial  velocities  are  set  outside  the  critical  range. 

<&pported  by  AF&ft  under  Grant  AFOSR-86-OI36A. 

I.  A.  Hirose  &  O.  ishihara,  Bull.  Am.  Phys.  Soc.  29,  1238 
(1984). 
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Qwjlliwif  Effect*  on  Turbulent  Bremsstr'ahlunt.*  O. 
Ishihara.  Texas  Tech  Univ..  and  A.  Hirose.  Unly/  of  Sas¬ 
katchewan— Th«  upconvereion  of  low  frequency  ion  acoustic 
turbulence  into  high  frequency  Langmuir  waves  b  re¬ 
examined  In  the  context  of  weak  turbulence  theory.  Con¬ 


trary  to  the  earOer  observations’*2  that  there  exists  turbu¬ 
lent  bremsstrahlung  of  Langmuir  waves  with  a  background 
of  Ion-acoustic  turbulence,  it  Is  shown  that  quasUinear 
effects  on  the  electron  dUtrlbutian  function  cause  stabiliza¬ 
tion  of  Langmuir  waves.  The  quesilinear  effects  are  found 
to  cancel  the  destabilizing  terms  originated  from  the  direct 
turbulent  corrections  to  the  Ibiear  response  taction.  To¬ 
gether  with  the  earlier  findings2*3  that  polarization  effects 
do  not  cause  destabilization  of  Langmuir  waves,  it  is  now 
concluded  that  there  should  be  no  upcoa version  of  ion 
acoustic  turbulence  into  Langmuir  turbulence,  which  agrees 
with  the  conclusion  derived  by  the  symmetry  argument/* 

•Supported  by  XFS5R  under  grant  AFOSR-46-0156A  (OX) 
and  NSERC,  Canada  {KMX 

1.  V.N.  Tsytovich,  L.  Stenflo  *  H.  Vlthefmsaon,  Physics 
Scripts  U,  251  (1975). 

2.  A.  Wrose,  Com.  Plasma  Phys.  Contr.  Fus.  L  117  (I9ML 

3.  D.F.  DufloU  &  O.  Pesme,  Phys.  FluMs  27.  2U  (I9ML 
.3.  Kuljpers  4  DA  Melrose,  Astrophys.  T  2H.  28  (1 98 5). 
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Stochastic  Haatlnt  of  Fluiu  In  Ion  Acoustic 

Turbulence.*  K.M.  TTJQI  and  0.  ISHIHAKA.  Tmi  Tech 
Unix-  Stochastic  ftocui  «£  plasma  tons  placed  In  the 
background  Ion  acoustic  Cuxbuloneo  if  ocudiod  analyti¬ 
cally  and  numarlcally.  Tbs  dispersion  relation  of  the 
loo  acoustic  w«v«  la  gluon  by  u  -  hC,.  uhoro  C,  Is  tha 
loo  acoustic  velocity.  tbs  randan  nature  of  tha  turbu- 
lonco  la  charactorlzod  by  tha  correlation  function  of 
alectrlc  fields.  ,We  nodal  tbo  turbulence  by  discrete 
Fourier  nodes  end'  tha  envelope  of  its  correlation  fune- 
don  ferns  a  Gaussian  spectrum.  -  By  controlling  the 
width  of  the  spectrun,  the  model  is  applicable  to  a 
single  node  as  a  limit  and  to  the  wide -band  turbulence. 
It  Is  shown  that  when  the  spectrum  width  bee  anas  wider 
tbs  basting  rata  of  plasma  particles  becomes  larger. 
Stochastic  heating  is  food  to  taka  place  for  loos  in 
resonance  and  also  not  la  resonance  with  Ion  acoustic 
waves. 


^Supported  by  APOSk  under  grant  APOSk-M-OlSdA. 
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ADI  2 

Plaema  Partlcla  Plffuslon  In  Prlft  Wove  Turbulence.  C. 
CIAMUSIU  sod  O.  1SH1HA8A,  This  Tech  P  *  Notions  of 
plasms  particles  have  been  studied  numerically  and 
analytically  when  the  particles  are  subjected  to  both  a 
background  magnetic  field  and  electric  fields  created  by 
drift  wave  turbulence.  For  these  studies,  test  parti¬ 
cles  are  placed  In  a  field  model  In  slab  geometry.  Both 
the  random  electric  fields  and  the  background  magnetic 
field,  with  Increasing  shear  la  the  x  direction,  are  set 
to  be  confined  in  the  y-s  plane.  Ifanerlcal  and  analyti¬ 
cal  results  indicate  that  particle  motions  increase 
along  Che  x  axis  as  the  magnitude  of  the  electric  fields 
Increases  and  as  the  density  of  the  mode  rational  sur¬ 
faces  Increases  along  the  x  axis.  «*Mre  tbs  mode  ratio¬ 
nal  surface  is  defined'  by  x  -  -n/u  with  n  and  m  being 
mode  number!  of  the  electric  fields. 


^Supported  by  AFOSK  under  grant  AFDSR-t(-01S(A. 


CONFERENCE  RECORD  -  ABSTRACTS 


1988  IEEE  INTERNATIONAL  CONFERENCE 
ON  PLASMA  SCIENCE 


June  6-8,  1988 
Seattle,  Washington 


4C6 

Stochastic  Diffusion  In  the  Absence  of  Resonances* 

0.  Ishlhara.  U.  Ho  and  K.M.  Yuen 
Department  of  Cleccrlcal  Engineering 
Texas  Tech  University.  Lubbock.  Texas 
and 

A.  Hlrose 

Departaent  of  Physics 

University  of  Saskatchewan.  Saskatoon,  Canada 


The  nonresonant  nature  of  charged  particle  notion 
in  the  electrostatic  plasma  turbulence  is  studied  mi* 
aerlcally  and  analytically  by  artificially  suppressing 
the  resonant  conditions  of  wave-particle  interaction. 
Nonresonant  structure  in  the  phase  space  is  realized 
in  the  test  particle  nunerlcsl  experiments  by  two 
different  methods.  One  approach  is  to  make  use  of 
discreteness  of  the  spectrum  of  plasma  turbulence. 
Test  particles  plsced  in  velocity  space  may  not  reach 
any  of  two  adjacent  modes  in  the  extremely  weak  turbu¬ 
lence.  This  situation  is  modeled  for  electrons  placed 
in  the  Langmuir  turbulence. *  Another  approach  Is  to 
place  test  particles  completely  out  of  the  resonant 
region  In  velocity  space.  This  situation  is  modeled 
for  ions  placed  in  the  Ion  acoustic  turbulence.2  Nu¬ 


merical  experiments  are  carried  out  for  test  particles 
placed  in  prescribed  random  electric  fields,  whose 
second  order  correlation  functions  are  characterized 
by  sinusoidal  structure  with  the  exponential  damping 
in  real  space,  described  by  the  summation  of  Fourier 
nodes  as 


E(x.t)  -  I  f(kj)cos(«(kj)t*kjX  ♦  6j). 


where  f(k«)  is  a  shaping  factor  with  Cauasian  form,  8j 
is  a  random  phase,  andwfkp  -  1  for  Langmuir  waves 
and  w(k<)  -  k«/(l*kj2)l/2  for  ion  acoustic  waves. 
Test  particles,  ranging  from  20  to  2000  in  number,  are 
followed  in  their  trajectories  in  the  discrete  fluctu¬ 
ation  spectrum  with  node  numbers  20  to  400.  Diffusion 
coefficients  are  measured  in  velocity  space  by  taking 
mean  square  variation  as  a  function  of  time  and  com¬ 
pared  with  the  analytic  predictions  evaluated  by 

0(v)  -^2/o"dt  /Jdx  «{x-vt)<£(x’.t*)E(x,*x.t'*t)>. 

where  e/m  is  a  charge /mass  ratio  and  set  to  be  1. 
Stochastic  diffusion  observed  numerically  both  In 
Langmuir  turbulence  and  in  ion  acoustic  turbulence  is 
in  good  agreement  with  the  theory  in  the  scale  of 
Kolmogorov  time  defined  by  tg  -  (k2  D)**/3.  In  the 
Langmuir  turbulence,  diffusion  is  observed  even  well 
below  the  critical  turbulence  level  given  by  the 
criterion  for  the  resonance  overlapping.  Our  observa¬ 
tion  suggests  the  possible  stochastic  acceleration  of 
charged  particles  in  the  nonresonant  nature  in  the 
plasma  turbulence. 

1.  0.  Ishlhara  and  A.  Hlrose,  Fhys.  Fluids  24.  2159 
(1985). 

2.  A.  Hlrose  and  0.  Ishlhara,  Bull.  Am.  Fhys.  Soc.  22. 
1238  (1984). 

*Vork  supported  by  AF0SR  under  Crant  AF0SR-86-01S6A 
(0.1 . .  V.H. .  K.Y.)  and  by  NSERC  of  Canada  (A.H.). 
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RESISTIVE  MHO  BALMXM1M  MOOC  IN  TORANAKS* 


A.  Hlroae 

Dept,  of  Physics.  Unlv.  of  Saskatchewan 
Sea Ratoon.  Saskatchewan.  S7N  0W0.  Canada 


0.  lahlhara 

Dept,  ef  Elect.  Enc. .  Texas  Tech  (tale. 
Lubbock.  TX  79409 


A  node  equation  far  the  realstive  MHO  ballooning 
■ode  la  lav  •  takaaak  nagnetlc  fsoaetry  baa  been 
derived  aad  analyzed  bath  analytically  and  numerically. 
Any  realatlve  aodes  In  a  hl(h  tanperatnra  tokanak 
rapilr*  ntranc  radial  locallzatloa  a boat  a  noda 
rational  onrface  la  order  to  satisfy  v(  » 

Hero,  v  la  tba  alectroa-loa  collision  frequency  and 
k.(r)  •  lur/L  la  tba  parallel  gradient  with  r  being 
the  radial  distance  froa  the  rational  so r face,  and  t 
the  sheer  length.  Strneg  radial  loca I last too  nakos  the 
nogaotlc  drift  frequency  operator 


*V.(ees*  - 


loins  »  ) 
k,  »rj 


secvlsr .  and  the  conventional  ordering 
w.  a,  »  «q 

totally  breaks  dew.  The  foraulatlon  In  the  present 
a  tody  is  free  froa  any  ordering  anong  the  frequencies 


v.  •  and  mq  and  thus  of fora  a  generalization  of  the 
conventional  node  equation  for  arbitrarily  strong 
radial  locallzatloa.  The  node  equation,  baaed  on 
electrostatic  Ion  response  end  folly  electromagnetic 
electron  response.  Is  given  In  Fourier  space  by 


where 


A  -  w 


D  «  «  -  ♦  Wj^IzMl**”1}  ♦  tjw  ♦  •ol«»  *  V 


0  •  Cw-#-wDtC«l](lt"1JwD#<«»*TCw-w0e(«»]Cw4w-il(l-r#| 

r#  -  e'bio(bj  .  b  -  (k,ptj*(i*«*)  .  T  «  \ni 


Mo  ens table  nodes  have  been  predicted  analytically  and 
found  nunerlcally  for  the  node  equation.  The  only 
bounded  solution  analytically  predicted  Is  the  stable 
Ion  dlanagnetlc  node.  The  denying  factor  Is  propor¬ 
tional  to  the  nagnetlc  sheer  and  *n.  the  toroldlclty-. 
This  has  been  fully  confirmed  by  ounerlcal  Integration 
with  para  asters  relevant  to  present-day  tokanaks. 


"Sponsored  by  RSERC.  Canada  (A.M.I  aad  AFOSR  (0.1.) 


Minutes  of  the  American  Physical  Society  Topical 
Conference  on  Plasma  Astrophysics 

Santa  Fe,  New  Mexico;  19-23  September  1988 


P212  Relativistic  Quasll Inear  Particle  Ac¬ 
celeration  In  Plasma  Turbulence.*  O.  ISHIHARA, 
Texas  Tech  P.  —Relativistic  quasllinear  kinetic 
theory  is  developed  to  describe  particle  accel¬ 
eration  in  the  presence  of  plasma  turbulence.  | 
Scattering  between  charged  particles  and  pies-  ; 
aons  is  considered.  It  is  assumed  that  each 
time  a  particle  makes  a  transition  the  energy 
it  gains  (or  loses)  comes  from  (or  goes  into) 
the  energy  of  turbulent  fluctuations.  The  theo¬ 
ry  includes  nonresonant  interactions  between 
charged  particles  and  plasaons  since  plasaons 
are  known  to  carry  not  only  electric  field  en¬ 
ergy  but  also  oscillating  energy  of  particle*.1 
The  formulation  takes  into  account  the  self 
fields  created  by  the  particle  itself,  which 
can  be  incorporated  as  the  effect  of  renormali¬ 
zation  of  the  particle  mass.2  The  derived  equa¬ 
tion  reduces  to  the  Balescu-Lenard  equation  in 
the  classical  limit.  A  possible  mechanism  of 
cosmic  ray  acceleration  is  discussed  based  on 
this  theory. 


♦Supported  by  OS  AFOSR  Grant  AFOSR-86-0156.. 
10.  Ishihara,  Pbys.  Rev.  A.  1£,  1219  (1967). 
2V.H.  Tsytovich,  Sov.  Phys.  JETP  51,  3  (1980) 
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Program  of  the  Thirtieth  Annual  Meeting 
of  the  Division  of  Plasma  Physics 
Hollywood,  Florida;  31  October-4  November  1988 


JV2I  Saturation  at  Won  linear  Drift  Ml 
*r  ToroldlcltV.*  o.  ISHXKASA,  Xnaui 
**  Pnivar-altv,  and  A.  BIROSE,  University  Of 
|M titrtunmn-  -Tohanaka  ara  known  to  bo  uneta- 
blo  with  respect  to  a  drift  mode  boooooo  of 
tboir  toroidicity.  The  turbulanoa  level  at 
saturation  of  tho  drift  soda  la  studiad  with 
tba  aaauaption  of  cold  ions  and  tba  alaotron 
tanparatura  gradient  as  a  principal  nonlinaar- 
lty.  A  positive  nonlinear  frequency  shift  is 
given  by  AQ  -  (v*2/12k2)  \4 I2#  where  fa  is 

the  alaotron  teaperature  gradient  relative  to 
a  density  gradient,  k  -  k|,a(k|  is  a  poloi- 
dal  azimuthal  va van  usher,  ,a  is  an  ion  Larmor 
radius  with  an  electron  tanparatura) ,  and  4  is 
the  peak  boundary  value  of  the  potential.  The 
nonlinearity  sanaitivaly  modifies  tba  radial 
profile  of  the  drift  mode,  which  in  tarn  in¬ 
duces  nonlinear  wave  damping.  The  resulting 
turbulence  level  at  saturation  falls  in  be¬ 
tween  the  strong  turbulence  limit  due  to  ion 
mode  coupling  and  the  weak  turbulanoa  limit 
due  to  electron  trapping. 

•Supported  by  HSERC,  Canada  (A.B.)  and  APOSR 
under  grant  S6-0156  (0.1.). 


6S 2  Tine-Dependent  Diffusion  Coefficient 
for  Hon-Markovlan  Process — In  Plassa  Turbu¬ 
lence.  *  H.  XIA  and  O.  XSHXHARA,  Texas  Tech 
University — Recent  study  of  resonance  broaden¬ 
ing  theory  showed  the  time  dependant  nature  of 
a  diffusion  coefficient  of  plasma  particles  in 
velocity  space,  2  where  a  particle  orbit  mod¬ 
ification  is  properly  treated.  Such  a  time 
dependence  of  diffusion  coefficient  is  studied 
based  on  the  generalised  Langevin  equation, 
where  a  frictional  force  exerted  by  the  effec¬ 
tive  collisions  due  to  plasma  turbulanoa  is 
included  with  its  retarded  effect.  This 
retarded  effect  is  rasponsible  for  the  non- 
Markovian  process  in  a  Gaussian  system.  It  is 
shown  that  a  derived  Fokker-Planek-type  equa¬ 
tion  has  a  diffusion  coefficient  which  depends 
on  time  and  reduces  to  the  conventional  quasi- 
1 inear  diffusion  coefficient  in  a  short  time 
limit  of  much  less  than  the  effective  colli¬ 
sion  time. 

•Supported  by  AFOSR  under  grant  SC-01SC. 

1.  0.  Xshihara  and  A.  Hirose,  Phys.  Fluids  28, 
2159  (19SS) . 

2.  A.  Salat,  Phys.  Fluids  31,  1499  (198S). 


Minutes  of  the  Seventh  Annual  Fall  Meeting 
of  the  Texas  Section 

Lubbock,  Texas;  4-5  November  1988 


turbulence  Is  studied 
end  analytically.  Turbaleuce  is  nodded  by 
random  electnc  fields  with  discrete  Fourier 
inodes  sad  the  conelatiou  function  of  the  fields 
is  assumed  to  be  Oeania.  Test  lots  sab  pieced 
in  s  background  of  ion  acoustic  fbdknoo  and 
the  statistics. of  particle  motions  ate  studied.  We 
have  observed  oat  stochastic  beating  can  take 
place  for  ions  both  in  resonance  mad  not  in 

reummee  with  inn  lemutie  waves.  1  1 


resonance  with  ion  nconsdc  waves.1 
particles  keep  absorbing  energy  from  the 
background  tuboleaee  while  the  noa-resooant 
particles  undergo  alternate  heating  and  cooling 
periods  with  n  net  heating  effect.  For  the*  ions' 
which  are  located  in  velocity  space  closer  to  the 
phase  velocity  of  the  ion  acoustic  waves  are 
observed  to  be  heated  more  effectively  than 
those  which  are  located  far  from  the  phase 

velocity  of  the  waves  .  The  numerical  results 
fairly  agree  with  the  analytical  prediction. 


*  Supported  by  AFOSR  under  grant  86-0156 
t.  KJM.  Yuen  and  O.  Ishihara,  ButL  Am.  Phys.  Soc. 
22,  1859  (1987). 


velocity  _  _  _ _  _ 

turbulence  was  presented  by  WdastocK1  and 
was  developed  ns  extended  resonance 
broadening  theories'1'3.  The  theories  were 
developed  oo  the  assumption  of  Gtasrina  and 
Markovian  stochastic  process.  Based  on  the 
Mori-Kubo  generalized  Lingerie  equation,  the 
exact  generalized  Fokker-Planck  conation 
corresponding  to  the  Gaussian  but  otherwise 
arbitrary  stochastic  force  is  studied.  It  agrees 
asymptotically  to  the  phenomenological 
Fcikker-Plaack  equation  for  systems  exhibiting 
short-time  tail  decay.  And  it  shows  n  velocity 
diffusion  coefficient  ns  n  function  of  time.  Tlmse 
theories  are  examined  by  the  numerical 
experieats  where  the  velocity  diffusion 
coefficients  are  measured  in  plasma  turbulence. 
*  Supported  by  AFOSR  under  nut  86-0IS6 

1.  J/Weiaoock.  Phys.  Fluids  12,  1045  (1969). 

2.  O.  Ishihara  and  A.  Hirose,  Phys.  Fluids  28, 
2159  (1985). 

3.  A  Salat,  Phys.  Fluids  31.  1499  (1988). 
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Particle  Diffusion  in  Turbulent 
Fields:  Transition  from  Quasilinear  to 
Nonlinear  Stage* 

Osaani  Xshihara 

Department  of  Electrical  Engineering 
Texas  Tech  University,  Lubbock,  Texas 

Charged  particles  in  randoei  electric 
fields  are  known  to  diffuse  over  phase 
velocity  resonant  layer  in  velocity  space 
(1,2],  while  charged  particles  placed  in  a 
sheared  magnetic  field  diffuse  over  the 
mode  rational  surfaces  due  to  SxB  random 
fluctuations  (3] .  He  study  in  detail  the 
transition  of  particle  diffusion  in 
turbulent  fields  from  quasilinear  (constant 
diffusion  rate)  to  nonlinear  regime  (tine 
dependent  rate)  by  test  particle  numerical 
experiments . 

We  consider  the  equation  of  motion: 

Atjdi  =£,  d£/dt  =  R(L0.  where  R(4.0  is  a  random 

electric  field,  ft-  C  R)  “  (*.  v,  qE/ta>)  for 
particle  motion  in  random  fields,  and 
(4.  4.  R)  *=  (y.  w  x  /L*.  cuEy/L«Bo)  for  guiding  center 
motion  in  a  sheared  magnetic  field,  where  Lj 
is  the  shear  length,  Hi  is  the  particle 
velocity  along  the  magnetic  field,  and  c  is 
the  speed  of  light.  Random  fields  are 
prescribed  to  have  Gaussian  structure  in 
Fourier  modes  with  the  spectrum  width  Sk. 
Trajectories  of  test  particles  are  followed 

and  statistics  of  (  variance  is  examined.  He 
define  the  diffusion  coefficient  by 

E^-idd^/dt  ,  where  the  time  slope  of  the  C 
variance  is  to  be  taken  at  t^t^te,  tK  is 
the  autocorrelation  time,  l/(£kt.),  and  tc  is 
the  time  at  which  particles  reach  the  edge 
of  the  resonance  layer.  He  have  observed 
that  particle  diffusion  rate  is  substan¬ 
tially  deviated  from  the  quasilinear  value 
anu  depends  on  time  when  the  amplitude  of 
random  field  becomes  larger.  Such  a 
transition  from  quasilinear  to  nonlinear 
diffusion  is  observed  without  particle  loss 
from  the  resonance  region.  Time  dependent 
diffusion  coefficient  may  be  explained  by- 
the  effect  of  retarded  friction  caused  by 
turbulent  fluctuations  (4]. 

1.  0.  Ishihara  and  A.  Hlrose,  Phys.  Fluids 
2A,  21S9  (1985). 

2.  A.  Salat,  Phys.  Fluids  21,  1499  (1988). 

3.  S.  P.  Hirshman  and  K.  Molvig,  Phys.  Rev. 
Lett.  12,  648  (1979). 

4.  H.  Xia  and  0.  Ishihara,  Bull.  Am.  Phys. 
Soc.  21,  2020  (1988). 

"Supported  by  the  AFOSR  under  Grant  86-0156. 
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Generalized  Langevin  Equation  for 
Plasma  Turbulence* 

Huajuan  Xia  and  Osamu  Ishihara 
Department  of  Electrical  Engineering 
Texas  Tech  University/  Lubbock,  Texas 

The  generalized  Langevin  equation, 
originally  derived  for  the  study  of 
Brownian  motion  by  Mori(l)  and  Kubo(2), 
has  been  applied  to  the  study  of  plasma 
turbulence.  The  evolution  equation  for  the 
particle  velocity  in  an  electrostatic 
plasma  turbulence  is  derived  by  the 
projection  operator  method.  The  time 
integral  appearing  in  the  equation 
indicates  that  the  wave-particle 
interaction  in  a  plasma  is  not  a  Markovian 
process.  Thus  the  wave-particle  interaction 
cannot  be  described  as  instantaneous 
collisions  between  particles  and  wave 
quanta  (quasi-particles)  as  is  usually 
done  in  a  conventional  weak  turbulence 
theory. 

He  consider  plasma  particles  placed  in 
an  electrostatic  plasma  turbulence 
characterized  by  random  electric  fields 
(Gaussian  process).  Particles  placed 
initially  at  (x,v)  move  in  the  random 
fields  to  (xt,  vt)  in  phase  space  at  time 
t.  Then  the  particle  velocity  can  be 
expressed  as 

vt  -  e  v 

where  -II  is  the  Liouville  operator  and 
is  given  by  -iJff  -  v  3/^x  «•  (e/*)E<x)?/?v. 
The  evolution  equation  (the  generalized 
Langevin  equation)  for  vt  can  be  derived  by 
the  method  of  the  projection  operator  as 

j  t 

=  -e’lQlE(x)  (x,v)  vt(x,v)  dt 

where  -iQ  -  -1(1-F)2J,  F  is  a  projection 
operator,  and  Yt_t  is  a  friction  function. 
The  time  integral  of  the  right  hand  side 
demonstrates  the  memory  effect  in  the 
wave-particle  interaction.  The  cumulant 
expansion  of  the  characteristic  function  of 
vt.  <exp(i£vt|>,  shows  the  time-dependent 
nature  of  the  velocity  variance.  Me  find 
that  the  memory  effect  or  the  influence  of 
the  past  on  particle  trajectories  in  phase 
space  limits  the  particle  diffusion  in 
velocity  space,  in  agreement  with  the 
earlier  observations  (3,4). 

1. H.  Mori,  Prog.Theor.Phys.32, 423  (1963). 

2. R.  Kubo,  Prog.Theor.Phys.22,  Part  1, 

255  (1966) . 

3.0.  Ishihara  and  A.  Hlrose,  Phys.  Fluids 
22,  2159  (1985). 

4. A.  Salat,  Phys.  Fluids  21,  1499  (1988). 

•Work  supported  by  AFOSR  under  Grant  86- 
0156. 
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PHOTON  ACCELERATION: 

THREE-WAVE  INTERACTION* 

_ O.  Ish&an  ! 

Department  of  Eloaricjl  Engineering 
Toms Tech  University,  Lnbbocfc.  Texas 

Reoendy  a  frequency  upshifoag  of  short  poises  of  User  I 

light  was  studied  by  Dawson  and  ocbcrsfl].  Pinions,  which 
propagate  in  a  plasma  with  effective  mass  of  txOpJc\  were 
shown  to  be  accelerated  (upshiftod  in  te)  by  propagating  down 
die  denafy  gradient  of  a  longibxfinal  plasma  wave. 

We  study  the  possibility  of  photon  aocdcradoa  from  die 
viewpoint  of  three-wave  interaction^  We  consider  that  a 

photon  with  momentum  fiq^  and  energy  fiOp.  will  combine 
a  plasmon  with  momentum  fiq^  and  energy  fwOpi  to  give  a 

photon  with  momentum  fMjpw  and  energy  The  waves 
will  interact  one  another  through  dm  wavoparticle  interaction  as  I 
shown  below: 

/  photon(f) 

I 

_ / 

/  / 

/  / 

photonG)  /  /  plasmon 

where  the  straight  line  indicates  a  plasma  particle  (charge  «, 
uw3s  tss  mid  P/.  The  iweMMUXs  vCmC*  uiCwicft  for 

the  pUsmon-panklecan  be  obtained  by  the  use  of  the  interaction 
Hamiltonian  between  particles  and  longitudinal  quasiparbcles 

Ht-|d3x'P+(x)e#(jt)'P(x)  . 

where  T(x)  is  a  wave  funaion  of  plasma  particles.  while  die 
interaction  vertex  function  for  the  photon -particle  can  be 
obtained  horn  the  interaction  Hanoi  Ionian  between  panicles  and 
transverae  qnasparades 

Hf  -mcf  **<*> <P  *  £*•>  *  A,(x) «P(x)  . 

where  A«  and  A|  are  the  zero-order  and  die  first-order  vector 
potential,  respectively .  Our  calculation  of  the  vertex  function  for 
this  three-wave  interaction  process  reveals  that  there  will  be  no 
transition  to  occur  if  the  plasma  is  cold  and  stationary.  We  have 
found  that  the  transition  probability  becomes  maximum  when 
die  plasma  panicles  move  in  die  direction  of  longitudinal  plasma 
oscillations  (pUsmoas)  along  with  the  direction  of  polamation 
of  vector  potential  A,  resulting  in  the  upshifting  of  frequency  of 
(he  outgoing  photon  fields. 

'Supported  by  AFOSR  under  Grant  86-0156. 

1 1 1  S.C  Wilks.  I.M.  Dawson,  W.8.  Mori.  T.  Katsouleas,  and 
M  E.  Jones.  Phys.  Rev.  Letters  62, 2600  (1989). 

|2|  E.G.  Harris,  in  Advances  in  Plasma  Physics,  ed.  by 
A.  Simon  and  W.B.  Thompson  (Interscicncc,  New  yotk, 
1969).  Vol.  3.  p.  157. 
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EFFECT  OF  TURBULENT  COLUSIONS  ON 
DIFFUSION  IN  STATIONARY  PLASMA  TURBULENCE* 


H.  XU  and  O.  Ishihara 
Department  of  Electrical  Engineering 
Texas  Tech  University,  Lubbock,  Texas 

Recently  the  velocity  diffusion  process  was  studied 
by  the  gcnetifiaedLangevio  equation  derived  by  the  projection 
operator  method!!].  The  further  study  shows  that  the  retarded 
frictional  function,  yft),  plays  an  important  role  in 
panicle  diffusion  in  the  velocity  space  in  stronger  tuibuleace  as 
much  as  the  resonance  broadening  effect.  The  retarded 
frictional  effect,  produced  by  the  effective  collisions  due  to  the 
plasma  imbulenoe,  is  described  as  j 

*«)— f - (Ex(x,v)JEt(x,v)) .  I 

nt^v.v) 


where  (.. . )  indicates  an  avenge  over  real  (x)  and  velocity  (v) 
spaces  and  random  fluctuation  field  ET(x,v)  is  assumed  to  be  a 
Gaussian,  but  noo- Markovian  and  non-wide-sense  stationary 
process.  The  velocity  diffusion  coefficient  in  a  stationary 
plasma  turbulence  in  the  absence  of  magnetic  field  is  given  as 


PfV.T„Tj. 


P(v.  th  Xi)  -^.(x.v),  E^fx.v)) 

«« X  Ck*-e»nXt|*xj)  (p-i  K4«,|-4«nl) 


with  <  >,  a  spatial  avenge,  U  the  resonance  broadening  term 
discussed  in  the  reference  2 .  Here 

zM-£fMi£k 

introduces  the  effect  of  turbulent  collisions  on  the  velocity  space 
diffusion  (tf»)  is  the  Fourier  transform  of  tft)).  The  diffusioa 
nte  becomes  time-dependent  because  of  the  non-Markovian  1 
nature  of  die  random  process  involved.  The  relation  between  , 
the  proposed  formulation  and  the  extended  resonance 
broadening  theory(2]  is  discussed.  We  also  carry  out  test  j 
panicle  numerical  experiment  for  Langmuir  turbulence  to  test  i 
the  theories.  In  a  stronger  turbulence  a  deviation  of  the 
diffusion  rate  from  the  one  predicted  by  both  the  quasilinear  and 
the  extended  resonance  theories  has  been  observed  and  is 
explained  qualitatively  by  the  present  formulation. 

*  Supponed  by  AFOSR  under  grant  86-0156. 

(1)  H.  Xia  and  O.  Ishihara.  Bull.  Am.  Phys.  Soc.  34. 1926 
(1989). 

(2)  O.  Ishihara  and  A.  Hirose,  Phys.  Fluids  2S,  2139  (1983). 


Program  of  the  Thirty-Second  Annual  Meeting 
of  the  Division  of  Plasma  Physics  .. 

Cincinnati,  Ohio;  12-16  November  1990 


ST 21  Photon  acceleration  hv  the  Plenmnn-Phnfrnn 
integration-*  o.  ishihara,  •»>*»<«  Tech  noiatc 
a  tty  -  Photon  acceleration  in  a  plasma  will  be 
observed  as  a  frequency  upahifting  in  the 
electromagnetic  wave  propagating  in  a  plasma. 
Me  study  photon  acceleration  from  the 
viewpoint  of  photon  emission  by  the 
interaction  between  plasmon  and  photon.  Three 
waves,  photon,  plasmon  and  photon,  interact 
one  another  through  the  interaction  between 
the  wave  and  plasma  particles.  The 
relativistic  interaction  Hamiltonian  between 
particles  of  species  a  and  photons  are  given 
by  H«-«'Id>*T+edTt  where  a  is  a  4X4  matrix 
expressed  in  terms  of  Pauli  matrices  and  a  is 
a  vector  potential,  while  the  interaction 
Hamiltonian  between  particles  and  plasmons  is 
H-e,ldJxT+4T,  where  4  is  a  scalar  potential. 
The  frequency  upshlfting  is  analyzed  in  tenss 
of  the  energy  gain  of  photon  as  a  result  of 
three  wave  interaction. 

♦Supported  by  AFOSR  under  grant  86-01S6. 
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2R1  Analysis  nf  the  Hon-Markovlan  Behavior  in 
Velocity  Diffusion  Process.*  H.  XIA  and  O. 
ISHIHARA,  Texas.  Tech  University  —  The  non- 
Markovian  diffusion  in  velocity  space  was 
observed  in  the  test particle  computer 
simulations  of  the  plasma  turbulence (1) .  The 
intrinsic  non-Markovian  effects  such  as 
resonance  broadening  and  retarded  turbulent 
collisions  were  found  to  be  responsible  for  such 
behavior [2] .  The  numerical  experiment  in  the 
moderately  strong  Langmuir  turbulence,  where  the 
resonance  region  is  bounded  by  the  lower 
boundary,  shows  that  particles  diffuse  toward 
the  upper  resonance  layer  after  they  hit  the 
lower  boundary.  It  is  found  that  the  non- 
Markovian  behavior  due  to  the  boundary  effect 
further  deviates  the  time-dependent  diffusion 
coefficient  predicted  by  the  intrinsic  non- 
Markovian  effects. 

•Supported  by  AFOSR  under  grant  86-0156. 

1.  O. Ishihara  and  A.Hirose,  Phys. Fluids  24,  2159 
<1985);  A. Salat,  Phys. Fluids  21,  1499  (1988). 

2.  H.Xia  and  O. Ishihara,  IEEE  Int.  Conf.  on 
Plasma  Science  (Oakland,  CA,  1990) . 


Abstract  Submitted  for  the 
18th  1991  IEEE  International  Conference 
on  Plasma  Science 


PHOTON  ACCELERATION 
BY  PLASMA  TURBULENCE* 


O.Ishihara 

IwnanniCOK  ol  UJOwUMhU  CX^EUKCuIIk 

Texas' Tech  University,  Lubbock,' Tacts 

Recently,  we  proposed  a  novel  method  of  photon 
icceleradont  in  which  clectnxntgDCtic  waves  injected  Into  a 
current  carrying  plasma  are  scattered  by  the  plasma  tnrtaleace 
and  dm  freqoency  of  the  electromagnetic  waves  are  upshlfted  by 
absorbing  energy  of  the  plasma  mrbulcttocj 

We  oandne  the  coaffltions  cf  efEecdvephoton  mwlpuiflii 
for  thf  above  meadooed  scheme  Cor  fats  pi»f»w 
panicles  ate  noa-reladvistk  or  relativistic.  We  Stan  with 
Hamiltonian  for  a  particle  of  species  s  in  an  electromagnetic 
Geld: 

l^"ilp;iAWP+**00 

i «  ca(p-^  A(x»  +  pmjC2 + e,6(x) 


where  A(x)  and  0(x)  are  the  vector  and  scalar  potentials  of  the 
Geld,  a  and  P  are  the  dimensionless  Dirac  4x4  matrices 
described  as  2  x  2  matrices.  The  three- wave  interaction  takes  : 
place  among  photon  (pj,  qO,  plasmon  (Qj,  qj)  and  photon 


(Qj,  qO-The  scattered  photoa  (Qj.  qj)  will  have  the  frequency 
of  Pj«Q|+ Pi  governed  by  the  conservation  law.  The  vertex 
function  for  die  photofrftasmonphoton  interaction  is  given  by* 


m  V  m - 

M-X  I 

•J  t*^W<^iQsl 

x  T-  f  *^*^«»«i^s(.. 

l  (Q,  -  v-qi)  (Qj  -  v-qj)  \Q,-vq, 


.1 


(Qj  -  v-qj)  (Q,  -  vq»)  \Qrv-qa  Prv-q, 


Qrv-q,/ 

)>. 


where  u,1#  and  u^  are  unit  polarization  vectors  of  photoa 
fields.  The  terms  with  n^  u^,  sre  resulted  from  relativistic 
effect  This  vertex  function  suggests  that  the  frequency 
uptUfkioc  becomes  suodiBom  amen  the  wave  vector  of  the 
emiue^gbotons  is  in  the  direction  of  electron  streams  in  the 
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Characteristics  of  Distribution  Function 
in  a  Time* Dependent  Velocity  Diffusion  Process* 

Husjuan  XU  and  Osama  Ishihara 

rv^nmimni  ofEbCtdai  Hnginwwiwg 

Texas  Tech  University,  Lubbock,  Texas 

The  time-dependency  of  die  particle  velocity  diffusion 
in  pi««m»  been  dbeuwed  by  nfftoH  d*wrwlcti 

modeis  and  shown  in  the  lest-paitide  computer  simulations  [1]- 
[4].  The  non-Maxkovian  effects  are  responsible  for  die  timo- 
depeiwent  nature  of  the  diffusion  coefficient  that  deviates  far 
from  the  quasiHncar  value  at  a  moderately  strong  turbulence 
level  in  Langmuir  turbulence  [5]. 

Under  the  assumption  of  Ocusslaa  turbulent  fields,  the 
oon-Markovian  effects  were  discussed  by  creating  the  panicle 
velocity  as  a  Wiener  process.  The  memory  property  appears  in 
the  spatial  Increment  of  particle  tndeokxy  as  well  as  In  me  tenn 
of  retarded  effective  turbulent  friction  la  the  generalized 
Langevin  equation  ($].  We  further  study  soch  a  noo-Marioovian 
behavior  in  vclodnrdigusionby  Intrododncthenon-Oanssian 
nature  in  nubulent  fields.  This  modification  leads  to  a  deviation 
of  probability  distribution  function  of  particle  velocity  from  the 
Gaussian  distribution.  The  effect  of  noo-Gaussian  turbulent 
fields  on  die  particle  diffusion  coefficient  Is  discussed  in  detafl. 

To  examine  the  long  time  behavior  of  the  diffusion 
process,  we  perform  test-particle  numerical  experiments  by  a 
tKwly  installed  MasPar  massively  parallel  computer  at  the  Texas 
Tech  University.  In  the  test-particle  simulation,  die  k-spectrum 
of  stationary  Langmuir  tmbulenoe  Is  modeled  by  two  different 
shapes;  one  is  Ga^enm-type,  and  the  other  Is  square-like. 
When  particles  hit  the  boundary  of  resonance  region,  a  fictitious 
non-Maikovian  effect  on  die  diffusion  coefficient  is  observed. 
By  varying  the  spectral  shape,  this  fictitious  noo-Markovian 
behavior  could  Se  suppressed.  The  results  of  test-particle 
numerical  experiments  am  presented  and  compared  with  the 
theory  based  on  die  assumption  of  noo-Oauwan  turbulent 
fields. 
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Nonresonant  wave-particle  interaction  in  sexniclassical  quasilinear  theory 

Osamu  Ishihara 

Department  of  Electrical  Engineering,  Texas  Tech  University,  P.O.  Box  4439,  Lubbock,  Texas  79409-4439 

(Received  20  August  1986) 

A  nonresonant  nature  of  wave-particle  interaction  is  cla rifted  from  the  viewpoint  of  quantum 
mechanics.  The  interaction  of  particles  and  quasi  particles  can  be  described  by  the  use  of  transition 
probability  which  is  found  to  have  both  resonant  and  nonresonant  contributions.  The  resonant  tran¬ 
sition  probability  is  known  as  Fermi’s  golden  rale,  which  is  now  supplemented  by  the  nonresonant 
contribution,  resulting  in  the  proper  conservation  of  energy  and  momentum  in  the  paxtide- 
quasi  particle  system. 
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DRIFT  ALFVEN  EIGENMODE  IN  TOKAMAKS 

A.  HIROSE 
Department  of  Physics, 

University  of  Saskatchewan, 

Saskatoon,  Saskatchewan, 

Canada 

O  ISHIHARA 

Department  of  Electrical  Engineering, 

Texas  Technical  University, 

Lubbock,  Texas, 

United  States  of  America 


A1STRACT.  A  drift  Alfda  skenmsils  having  a  phaaevalocltyalaBlQcaatly  larger  tha  aWaw.-nf— 

velocity  h*own  to  exist  In  the  tnlramatr  magnetic  gsumstt|.  For  typlesl  ohmlcal  heated  Soksmsk  discharges  the 
sppraxlasla  flap  infos  nhdoa  h  warScj.,.  Rich  a  Zmt  drift. amds  has  racentty  bean  nbssraad  la  the  TEXT 
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SEARCH  FOR 

ION  TEMPERATURE  GRADIENT  DRIVEN 
ELECTROSTATIC  HYDRODYNAMIC  INSTABILITY 
IN  TOKAMAKS 

A.  HIROSE 
Department  of  Physics, 

University  of  Saskatchewan, 

Saskatoon,  Saskatchewan, 

Canada 

O.  ISHIHARA 

Department  of  Electrical  Engineering, 

Texas  Tech  University, 

Lubbock,  Texas, 

United  States  of  America 


ABSTRACT.  The  Ion  temperature  gradient  driven  electrostatic  instability  in  tokamak  magnetic  geometry  has 
been  investigated  numerically  in  terms  of  a  generalized  ion  density  moment  For  typical  tokamak  parameters 
(ratio  between  the  density  gradient  scale  length  and  the  major  radius,  e,  *>  L„/R,  about  0  (0.1)  and  T,e;Tt)  the 
search  for  rapidly  growing  instabilities  of  a  hydrodynamic  nature  has  been  unsuccessful,  even  with  the  ion  tempera¬ 
ture  gradient  relative  to  the  density  gradient,  =  d(InT,)/d(lnn0),  as  large  as  six. 
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Resistive  ballooning  mode  in  tokamaks 

A.  HotOSE  andT.  L.  Kroeker1 

Department  ef  Physics,  University  of  Saskatchewan,  Saskatoon,  SasL,  Canada  S7N  0W0 

AND 

O.  ISHIHAKA 

Department  of  Electrical  Engineering,  Texas  Tech  University,  Lubbock,  TX  79409.  USA. 

Received  May  3, 1988* 

The  resistive  ballooning  mode  in  tokamak  geometry  has  been  rc-invcstigaacd  without  imposing  any  ordering  among  cfarac- 
tesMc  fieqwades,  |**j  (growth  ate),  at*  (diamagnetic  frequency),  and  **d(r)  (magnetiodrift  frequency  fat  die  ballooning 

apace)  *•*««■■  -Ti  afcal - h htashMe — **‘ni  fin  a itriilr  log  ifiniimgnrtir  mnrtr  (analytically 

profioed)  hat  been  found. 
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RESISTIVE  MHD  BALLOONING  MODE  IN  TOKAMAKS 

A.  HIROSE 
Department  of  Physics, 

University  of  Saskatchewan, 
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Canada 

O.  ISHIHARA 

Department  of  Electrical  Engineering, 

Texas  Tech  Univerrity, 

Lubbock,  Texas, 

United  Slates  of  America 


ABSTRACT.  A  Bode  eqvadoa  far  the  resistive  MHD  taDoadif  mode  in  low  beat  tokxmak  magnetic  geometry 
has  been  derived  and  analysed  »Mwa  imposing  any  ordering  among  (he  iraqucada  a,  «,  (diamagnetic  froqoeacy) 
and  ug  (magnetic  drift  freqoency  operator)-  Only  a  stable,  ion  AMipmii'  mode  has  been  found.  The  verisrive 
ballooning  mode  is  «tro«uty  stahilimil  by  magnetic  risear  and  toroidal  effects. 
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Resonance  broadening  in  drift  wave  turbulence 


O.  Ishihara  and  C.  Grabowskr*’ 

Department  of  Electrical  Engineering.  Texas  Tech  University.  Lubbock,  Texas  79409 
A.  Hirose 

Department  of  Physics,  University  pf Saskatchewan.  Saskatoon.  Saskatchewan.  S7N  OtVO,  Canada 
(Received  30  June  1989;  accepted  4  October  1989) 

Turbulent  diffusion  of  electrons  in  drift  modes  in  a  sheared  magnetic  field  is  studied  by  using  a 
test  particle  numerical  experiment  Electrons  diffuse  across  the  magnetic  field  over  mode 
rational  surfaces,  where  elections  interact  with  waves  in  resonance.  A  spatial  diffusion 
coefficient  which  describes  resonance  broadening,  is  found  to  be  time  dependent  and  departs 
from  quasilinear  predictions  in  stronger  turbulenoe  even  well  before  the  time  when  particles  hit 
the  resonance  boundary  of  the  rational  surfaces. 
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Photon-Plasmon-Photon  Interaction 

Osamu  Ishihara 

Department  of  Electrical  Engineering  and  Department  of  Physics 
Texas  Tech  University,  Lubbock.  Texas 79409-3102 


ABSTRACT 

A  theory  of  photoo-plasmon-photon  three-wave  interaction  Is  presented.  An  electromagnetic 
wave  launched  into  the  plasma  turbulenoe  interacts  each  other  through  die  wave-panicle  interaction 
resulting  in  the  emission  of  an  electromagnetic  wave  with  frequency  upshifted.  The  theory  is 
based  on  the  Dirac  relativistic  wave  equation,  where  the  interaction  of  plasma  particles  with 
electromagnetic  waves  are  described.  It  is  shown  that  die  frequency  upshifting  of  the 
electromagnetic  wave  can  take  place  in  a  current  carrying  plasma. 


PACS  numbers:  52.40.Db,  52.35.Ra,  52.40.Nk 


